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ABSTRACT 

Using the POLISH instrument, I am unable to reproduce the large-amphtude 
polarimetric observations of Berdyugina et ah (2008) to the > 99.99% 
confidence level. I observe no significant polarimetric variability in the HD 
189733 system, and the upper limit to variability from the exoplanet is 
AP < 7.9 X 10^^ with 99% confidence in the 400 nm to 675 nm wavelength 
range. Berdyugina et al. (2008) report polarized, scattered light from the 
atmosphere of the HD 189733b hot Jupiter with an amplitude of two parts 
in 10^. Such a large amplitude is over an order of magnitude larger than 
expected given a geometric albedo similar to other hot Jup iters. However, my 
non-detection of polarimetric variability phase-locked to the orbital period of the 
exoplanet, and the lack of any significant variability, shows that the polarimetric 
modulation reported by Berdyugina et al. (2008) cannot be due to the exoplanet. 



Subject headings: instrumentation: polarimeters — planetary systems — polar- 
ization — stars: individual (HD 189733) — techniques: polarimetric 



1. Introduction 

Radial velocity surveys uncover large populations of exoplanets that allow models of 
planet formation and migration to be constrained. However, in order to study individual ex- 
oplanets, it is important to directly detect both their scattered optical fiux and their thermal 
emission. While the Spitzer Space Telescope has enabled exoplanetary thermal emission to 
be observed (Deming et al. 2005; Harrington et al. 2006, 2007; Knutson et al. 2007, 2009a, 
2009b), scattered fiux from only one exoplanet has been conclusively observed (Kalas et al. 
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2008). This is because the contrast ratio between star and exoplanet is at least an order 
of magnitude larger in the optical than in the infrared. Thermal emission measurements 
allow exoplanetary temperature maps to be made (Knutson et al. 2007), which constrain 
models regarding redistribution of stellar insolation by exoplanetary winds. In addition, 
detection of infrared emission allows molecules such as water vapor (Tinetti et al. 2007) 
and methane (Swain et al. 2008) to be identified in exoplanetary atmospheres. Atomic 
species can be identified in these atmospheres by transmission spectroscopy in the optical, 
and it can also be used to identify the presence of cloud layers (Charbonneau et al. 2002; 
Vidal-Madjar et al. 2003; Barman 2007; Pont et al. 2008; Redfield et al. 2008). If the 
exoplanetary radius can be estimated, detection of light scattered by an exoplanet allows 
its geometric albedo to be determined. Geometric albedo is a measure of the scattering in 
the atmosphere of the exoplanet, which gives information about atmospheric cloud structure. 

However, most thermal emission and transmission spectroscopy measurements of exo- 
planets are from transiting systems, where the orbital plane is seen edge-on. This is because 
thermal emission may be identified during secondary eclipse, when the exoplanet is occulted 
by the star, and transmission spectroscopy requires a primary eclipse by definition. Tran- 
sits allow baseline stellar emission (both in the optical and infrared) to be subtracted from 
the combined star/exoplanet signal, which greatly improves the signal to noise ratio for di- 
rect detection. Unfortunately, transiting systems only comprise about 10% of massive, short 
period exoplanets (so-called "hot Jupiters" ) , so high signal to noise observations of the atmo- 
spheres of most known exoplanets are not accessible with these techniques. While advanced 
imaging instruments, such as the Gemini Planet Imager (Macintosh et al. 2006; Graham et 
al. 2007b), have the potential to observe thermal emission from long-period exoplanets, the 
majority of known exoplanets orbit too closely to their host stars in order to be accessible 
to imaging. 

In addition to observing scattered and emitted fiux from exoplanets, determining ac- 
curate exoplanetary masses is needed to characterize individual exoplanets. However, the 
dominant exoplanet-finding technique, radial velocity, is insensitive to stellar reflex motion 
in the plane of the sky. Therefore, measured mass m is only a lower hmit to the true mass 
M because m = Msini. While one may infer the most probable exoplanetary mass by 
assuming an isotropic distribution of orbital inclination i, an observational technique that 
constrains inclination is desired. 

Again, transiting systems are a boon to the characterization of individual exoplanets. 
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because the shape of the transit hghtcurve is indicative of both orbital inchnation and exo- 
planetary radius. Inchnation estimates from transit observations can be coupled with radial 
velocity data to derive accurate masses. Indeed, the mass of the transiting hot Jupiter HD 
209458b has been measured to within five Earth masses (Torres et al. 2008 analysis of Brown 
et al. 2001; Naef et al. 2004; Butler et al. 2006). The hot Neptune GJ 436b, with an M 
dwarf host, has an uncertainty of less than one Earth mass (Torres et al. 2008 analysis of 
Deming et al. 2007; Gillon et al. 2007a, 2007b; Maness et al. 2007). Knowledge of accurate 
exoplanetary radii enables bulk exoplanetary density to be determined. For example, the 
relatively large mass with respect to radius of the transiting exoplanet HD 149026b indicates 
the presence of a large, rocky core, which bolsters support for its formation via core accre- 
tion (Sato et al. 2005). In order to determine accurate exoplanetary mass for non-transiting 
exoplanets, however, new observational techniques are necessary. 



2. Exoplanetary Polarimetry 

Polarimetric observations have the potential to determine orbital inclination, and there- 
fore accurate mass, for exoplanets by isolating their scattered, optical flux. Indeed, po- 
larimetry has been utilized to study dust grains in debris disks because of its sensitivity to 
scattered hght (Meyer et al. 2002; Jensen et al. 2004; Monin et al. 2006; Perrin et al. 2006; 
Graham et al. 2007a; Beckford et al. 2008). The extra information content offered by po- 
larimetry over photometry suggests that exoplanetary scattered flux can be detected even in 
non-transiting systems. That is, photometric detection of exoplanetary flux benefits greatly 
from chopping of exoplanetary flux between in-transit and out of transit observations due to 
the scalar nature of photometric observations. However, the vector nature of polarimetric 
observations may enable direct detection of exoplanetary flux even when system intensity is 
constant. 

Polarization of exoplanets arises by scattering of incident starlight by gas molecules, 
aerosols, and dust grains in the atmosphere of the exoplanet (Seager et al. 2000). For a 
face-on, circular orbit (Figured^), the exoplanet is always seen at quarter phase and always 
has half its disk illuminated. For a featureless exoplanetary atmosphere, the intensity of light 
scattered by the exoplanet is constant throughout the orbit and the degree of polarization is 
also constant. However, the position angle of polarization rotates through 360° each orbit, 
because the scattering plane rotates as the exoplanet progresses in its orbit. 
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Fig. 1. — Schematic orbital modulation of polarization for an exoplanet with (a) face-on and 
(b) edge-on geometries. The amount of light scattered by the exoplanet is represented by 
the white, illuminated portion of the exoplanet, and the degree of polarization is strongest 
when the exoplanet is near quarter phase. The position angle of net polarization is given by 
the orientation of the black lines. 

In contrast, an edge-on viewing geometry generates large, periodic variability in the de- 
gree of polarization because the amount of scattered light is variable (Figure [T)d). However, 
the scattering plane is always nearly coplanar with the orbital plane, so the position angle 
of net polarization does not vary significantly throughout the orbit. It should be noted that 
the polarimetric signature of an exoplanetary transit is discussed in Carciofi & Magalhaes 
(2005) and in section 4.2. In general, an exoplanet exhibits variability in the polarization 
vector that is indicative of orbital inclination, and the models of Seager et al. (2000) and 
Stam et al. (2004) demonstrate this effect. 

The simplest discussion of the modulation in system polarization due to an exoplanet 
exists for a Lambertian phase function and a polarization dependence on scattering angle 
similar to Rayleigh scattering. This is given by 
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P{ct>) = eF{ct>)P,{ct>). (1) 

Here, is orbital phase pinned to stellar radial velocity phase (0 = represents superior 
conjunction of the exoplanet), e = p^Rp/af is the fraction of stellar flux scattered by an 
exoplanet with radius Rp and semimajor axis a (assuming a circular orbit), p is its geometric 
albedo (the fraction of exoplanetary scattered flux at full phase compared to that scattered 
by a Lambertian disk), -F(0) is the phase function (the fraction of intercepted flux scattered 
toward the observer), and -Po(0) is the polarization of that scattered flux. The phase function 
is given in terms of the phase angle a, which is the angle between the host star and observer 
as seen from the exoplanet. Russell (1916) gives the Lambertian phase function as the 
following: 



Fia) 



sin a + [tx — a) cos a 



TT 



(2) 



Given cos a = sin i cos 0, -F(0) can be determined. It is useful to decompose -P(0) in terms 
of its normalized Stokes parameters, where -P(0) = a/ Q^(0) + ^^(0) and 

Q'(0) = eF(0) (sin^ - cos^ cos^ i) (3a) 
U\<P) = eF{<p) sin 20 cos i (3b) 



after Shakhovskoi (1965). The primed Stokes parameters are measured in the orbital frame 
and +Q' is defined to be in the direction of the orbital angular momentum vector. Con- 
verting polarization from the orbital frame to celestial coordinates requires knowledge of the 
longitude of the ascending node where = VL + 270° (Figure [2]): 

Q(0) = g'(0) cos 29 - U'{(j)) sin 26 (4a) 
U{(f)) = g'(0) sin 26 + t/'(0) cos 26. (4b) 



While terrestrial planets are expected to follow AP = P(0)max — -P(0)min oc p because of 
their solid scattering surfaces, multiple scattering reduces the polarization of light scattered 
by gas giant atmospheres with high albedos (Seager et al. 2000; Stam et al. 2004; Lucas et 
al. 2009). Thus, the form of polarization as a function of orbital phase from Equations 1 
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Fig. 2. — Schematic polarization due to an exoplanet at quarter phase and to an equatorial 
starspot on the stellar limb (section 4.2). Celestial North is in the direction of +Q, and East 
is in the direction of —Q. Sizes of the star, exoplanet, and orbit are to scale for the HD 
189733 system. The value of = 16° is taken from Berdyugina et al. (2008). 

through 4 is expected to approximate that of all exoplanets, but the amplitude of the signal 
will be strongly overestimated for gas giants. 

From the Monte Carlo scattering simulations of Lucas et al. (2009), the polarization 
amplitude of a gas giant exoplanet is expected to be 

AP..,.. = (-2.10/ + 2.91rt (^) X 10- (5a) 

AP^ij = 1.43(APi.s„.). (5b) 



Equation [5a| is the polarization amplitude for an edge-on orientation, while Equation [5b] 
relates to the general case where orbital inclination is unknown. The quadratic fit to AP as 
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a function of p reproduces the <Pmv> values in Table 6 of Lucas et al. (2009) to within 2%. 

The strongest polarimetric modulation from a gas giant exoplanet will occur for a Lam- 
bertian phase function with a polarization dependence on scattering angle similar to Rayleigh 
scattering. Here, p = 2/3 at all wavelengths. Since {Rp/af = (3.081 ± 0.043) x 10"^ for HD 
189733b (Torres et al. 2008 analysis of Knutson et al. 2007), Equation [5al shows that the am- 
plitude of polarimetric modulation due to the exoplanet cannot be larger than AP ^ 2x 10~^ 
for scattered light at all wavelengths. Upper limits to geometric albedo from reflected light 
campaigns exist for the exoplanets around r Boo, HD 75289, and HD 209458 from photom- 
etry (Leigh et al. 2003a, 2003b; Rowe et al. 2008b). Lucas et al. (2009) estimate an upper 
limit for r Boo b using PlanetPol, a polarimeter capable of part per million precision on 
bright stars (Hough et al. 2006). These albedo upper limits, along with the bandpass in 
which they were determined, are shown in Table [H Assuming a mean albedo upper limit 
of p < 0.22, the polarimetric modulation from HD 189733b is APcxp ^ 1 ^ 10"^. Since the 
albedo estimates in Table [1] range from 400 to 920 nm, this expected upper limit to HD 
189733b modulation is for scattered light at all optical wavelengths. 



3. Observations and Data Reduction 

The POLISH instrument (POLarimeter for Inclination Studies of High mass x-ray bi- 
naries/Hot jupiters) is a visible light polarimeter commissioned at the Cassegrain focus of 
the Hale 5-m telescope at Palomar Observatory, California. This instrument is described 
in Wiktorowicz & Matthews (2008, hereafter WM08). It utilizes a photoelastic modulator 
(PEM), a Wollaston prism, lock-in amplifiers, and digital voltmeters to modulate and de- 
tect incident, polarized light at 100 kHz. The Wollaston prism feeds a pair of avalanche 
photodiodes (APDs) or photomultiplier tubes (PMTs), depending on stellar intensity. The 
bandpass of the instrument is limited by the detectors and atmospheric transmission; the 
lack of spectral filters increases throughput of the instrument and allows for high precision 
observations. The APD bandpass ranges from about 400 to 850 nm, while the PMT band- 
pass ranges from roughly 400 to 675 nm. Atmospheric transmission sets the lower limit 
for spectral coverage for both detector types. The above components and large telescope 
aperture contribute to the high signal-to-noise observations with the instrument, where part 
per million precision is achieved on bright stars (WM08). On-source guiding is accomplished 
by use of a beamsplitter, which allows about 5% of the flux to be sent to a CCD camera. 
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The voltage output from the detectors consists of a roughly sinusoidal waveform with 
a frequency of 100 kHz and a DC offset. Polarization of the input beam is proportional 
to the ratio of the AC amplitude (polarized flux) to the DC offset (total flux). Thus, at- 
mospheric effects that operate on timescales longer than 10 /is arc effectively removed. For 
example, non-birefringent cirrus clouds passing through the telescope beam will reduce both 
polarized and unpolarized flux equally: the ratio of AC to DC signals will remain constant. 
Additionally, polarization imparted from such forward scattering by cloud particles will be 
zero. Therefore, high precision, integrated hght, polarimetric observations with POLISH 
are photon shot noise hmited (WM08), eliminating the benefit from a similar, space-based 
polarimeter. 

The amplitude of the 100 kHz AC signal is sampled at 8 Hz by the lock-in amphfiers, 
and the DC offset is sampled at 6 Hz by the voltmeters. The fact that the sampling rates 
are different is not important because polarization is related to time-averaged AC and DC 
values. Each on-source measurement consists of one, 30 second integration during which the 
AC amplitude and DC offset are sampled. AC and DC signals are then sky subtracted by 
chopping the secondary mirror to a sky field and initiating another 30 second integration. 
It should be noted that no less than 30 seconds elapse before the next chop throw. Since 
sky background is low in both the AC amplitude and DC offset, asymmetric, 2:1 source/sky 
chops are used. That is, observing "triplets" consist of repeated sequences of source, sky, 
and source integrations every 90 seconds. 

The Cassegrain ring is rotated by 45° roughly every 10 minutes, which causes a different 
linear Stokes parameter to be sampled. This ensures that all ±Q and ±U Stokes parameters 
are sampled at least once for each star to minimize systematic effects. Note that each star 
is observed for about one hour per night. Nightly telescope polarization of 3.0 x 10^ is then 
subtracted. Telescope polarization is the mean of the nightly measurements of the bright 
(y ~ 3), unpolarized star HR 5854 weighted by the DC level of each measurement. DC level 
is proportional to the number of detected photons, and weighting by this quantity ensures 
that each photon is weighted equally. This is important in partly cloudy conditions, because 
observed stellar intensity varies throughout the night. 

The equatorial mount of the Hale 5-m inhibits traditional telescope polarization mea- 
surement, which involves allowing the field to rotate and determining the center of the {Q, U) 
locus for weakly polarized targets. However, Hough et al. (2006) and Lucas et al. (2009) 
performed this analysis and report a polarization of four parts per million for HR 5854. Lucas 
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et al. (2009) report that the polarization of this star varies by seven parts per milhon over 
one year, but this hes at the < Sa level of significance. Therefore, I assume that this star is 
indeed unpolarized, and I assume that it is non- variable over this seven-night observing run. 
Since telescope polarization is a function of wavelength, it will be different for stars with 
different spectral types. In addition, I quickly measure telescope polarization to the part per 
million level with APDs on a bright, unpolarized star, even though the science targets in 
this work are about five magnitudes fainter and are observed with the different bandpass of 
the PMTs. Therefore, an error of order one part in 10^ or less in the absolute calibration of 
telescope polarization may be present. However, such an offset is assumed to be non- variable 
and below the measurement noise for the targets in this work. Systematic effects have been 
reduced to the part-per-million level for weakly polarized sources and to less than 1% of the 
measured polarization for significantly polarized sources (WM08). 

Nightly mean and run-averaged polarizations for each target are determined by tak- 
ing the mean of the polarization measurements weighted by their DC levels. Uncertainty 
in polarization is the square root of the weighted variance of the polarizations from these 
measurements divided by the square root of the number of measurements. Each polariza- 
tion measurement is photon shot noise limited, and uncertainty in run-averaged polarization 
scales with the square root of stellar polarization. The polarization noise floor of the instru- 
ment is about two parts per million (WM08). I calibrate absolute polarization measured 
with POLISH against the catalog of Heiles (2000). However, since my observed position 
angles of net polarization are generally consistent with those in the literature for strongly 
polarized targets, I do not calibrate position angle. It should be stressed that the unfiltered 
polarimetric observations of science targets in this work are taken over the wavelength range 
of roughly 400 to 675 nm. 

The stars observed are listed in Table [2j Polarimetric data are obtained from the cat- 
alogs of Heiles (2000), from Berdyugina et al. (2008, hereafter B08), and from WM08. All 
other non-polarimetric data are taken from the SIMBAD database. The polarization and 
position angle values in parentheses represent the standard error of the mean, which assumes 
nonvariable target polarization. However, Cygnus X-1 is known to be variable with a po- 
larimetric amplitude of order AP ^ 0.1% (Nolt et al. 1975; Kemp 1980; Dolan & Tapia 
1989). The amplitude of variability of Cygnus X-1 is also variable and due to both stable 
and stochastic effects (Dolan & Tapia 1989, 1992). This variability is roughly phase-locked 
to the orbital period of the binary, and Cygnus X-1 is included in this work as a variable 
control source. 



-10- 



Table 1. Upper Limits to Exoplanetary Albedos 



Target 


P 


Confidence 






Method 


Ref 








(nm) 


(nm) 






HD 75289 


< 0.12 


99.9% 


406.5 


522.0 


Photometry 


1 


HD 209458 


< 0.17 


3(7 


400 


700 


Photometry 


2 


r Boo 


< 0.37 


4a 


590 


920 


Polarimetry 


3 


r Boo 


< 0.39 


99.9% 


407.4 


649.0 


Photometry 


4 



References. — (1) Leigh et al. (2003b), (2) Rowe et al. (2008b), (3) 
Lucas et al. (2009), (4) Leigh et al. (2003a). 



Table 2. Observed Stars 



Name 


Alt. Name 


RA 


Dec 


P 


e 


Ref 


V 


Type 






(J2000) 


(J2000) 


(%) 


n 








HR 5854 


a Scr 


15h44m;^gSQy 


+06?25'32('3 


0.00020(22) 




1 


2.64 


K2inb 


HD 149026=' 


SAO 65349 


16'' 30™ 29! 62 


+38?20'50'.'3 


7 


? 




8.16 


GOIV 


HD 175541*' 


GJ 736 


18''55™40':88 


+04?15'55'.'2 


7 


? 




8.03 


G8V 


Cygnus X-l<= 


SAO 69181 


19''58™21?68 


+35?12'05'.'8 


4.98(18) 


135.0(1.0) 


2 


8.95 


09.7Iab 


HD 189733'' 


V452 Vul 


20''00™43?71 


+22 ? 42' 39 '.'1 


0.0212(30) 


99.6(4.1) 


3 


7.68 


K1.5V 


HR 8974^' 


7 Cep 


23''39™20!85 


+77?37'56'.'2 


0.00047(35) 




1 


3.23 


KlIV 



'"Transiting, short period exoplanet host 
''Non-transiting, long period exoplanet host 
"High mass X-ray binary 



References. — (1) WM08, (2) Heiles (2000), (3) B08 

Note. — WM08 polarizations are measured in the wavelength range of about 400 to 850 nm, Heiles (2000) 
polarization data are quoted for V band, and B08 polarization is measured in B band. 
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The long period exoplanet hosts HD 175541 and HR 8974 are included as non-variable 
control sources, as the exoplanets should only be detectable over an entire orbit and at 
the level of APexp < 2 x 10"^ and APexp < 5 x 10"^ respectively (Equation [5b]). The 
transiting hot Saturn around HD 149026 is expected to generate a polarimetric signal of 
APexp ^ 2 X 10~^ over its orbit (Equation [5a]). Orbital parameters and expected polariza- 
tions of these systems are given in Table [31 These amplitudes assume a geometric albedo 
upper limit of p < 0.22 (the mean of the upper limits in Table [1]) and a radius of P = 1.2 
Rj for the non-transiting exoplanets. Since these upper limits to polarimetric modulation 
range from 400 to 920 nm, the expected upper limits are for scattered light at all optical 
wavelengths. No polarimetric data for HD 149026 or HD 175541 were found in the literature. 



4. Results and Discussion 

4.1. Unbiased Significance of Variability 

Nightly mean and run-averaged polarizations for the observed targets are listed in Table 
[H after telescope polarization (hereafter "TP" ) is subtracted. Note that run-averaged po- 
larization assumes non-variability of the source. This is apparent in the formal uncertainty 
in polarization for Cygnus X-1, ap ~ 4 x 10~^, because this source is known to be variable 
on the order of AP ^0.1%. I refer the reader to Table 4 of WM08 for nightly polarization 
measurements of the unpolarized standard star HR 8974. The standard deviation of mean 
nightly polarizations is given as APobs in Table [31 These represent measurement precision 
coupled with intrinsic source variability. 

To determine the unbiased significance of night-to-night variability of each target, I 
perform a Kolmogorov-Smirnov (K-S) test on the measurement distribution for each pair 
of nights. This test estimates the significance of variability without prior knowledge of its 
physical cause. This is because analysis of data from all pairs of nights samples all temporal 
frequencies available as opposed to focusing on the frequency of an expected signal. For 
instance. Equations [H through 4 indicate that polarization due to stellar flux scattered off an 
exoplanetary atmosphere will be strongest at quarter phases and will be zero at conjunctions. 
Thus, the unbiased K-S test will underestimate the significance of variability caused by this 
process. However, this test is important in separating stochastic variability of the host star 
from exoplanetary modulation. 
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Table 3. Expected and Observed Variability 



R a z Ref APexp 



(^j) (Au) n 



Target 



HR 8974 
HD 175541 
HD 149026 
HD 189733 
Cygnus X-1 



«i 1.2 
K 1.2 

0.654(;60) 

1.138(27) 



2.044(57) 
1.03(-) 

0.04313(1^^) 
0.03099(1^11) 
0.195(42) 



? 1 

? 2 

ntt-'o) 3 

85.58(6) 4 

48(7) 5 



< 5 X 10-9 

< 2 X 10-8 

< 2 X 10-6 

< 1 X 10-5 
« 10-2 



APobs 



9.0 X 10-6 
2.6 X 10-5 
1.6 X 10-5 

2.1 X 10-5 
5.1 X 10-"^ 



References. — (1) Neuhauser et al. (2007); (2) Johnson et al. (2007); (3) Torres 
et al. (2008) analysis of Sato et al. (2005), Butler et al. (2006), Charbonneau et al. 
(2006, Winn et al. (2008); (4) Torres et al. (2008) analysis of Bouchy et al. (2005), 
Knutson et al. (2007); (5) lorio (2008) analysis of Gies et al. (2003), Shaposhnikov & 
Titarchuk (2007). 

Note. — Expected variability for all targets is in broadband BVRI light. The 
observed upper limit to variability of HR 8974 is in the wavelength range of about 
400 nm to roughly 850 nm, and upper limits for the rest of the sample are measured 
between 400 and 675 nm. References are for orbital data. 
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Table 4. Observed Polarization 
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2008 Jun 08 


HD 189733 


-0.0174(42) 


+0.0190(31) 


0.0258(37) 


66.3(4.2) 


2008 Jun 09 




-0.0255(12) 


+0.0172(17) 


0.0307(14) 


73.0(1.4) 


2008 Jun 10 




-0.0206(15) 


+0.0125(18) 


0.0241(16) 


74.4(2.1) 


2008 Jun 11 




-0.0209(17) 


+0.0176(17) 


0.0274(17) 


70.0(1.8) 


2008 Jun 12 












2008 Jun 13 




-0.0242(12) 


+0.0163(16) 


0.0292(14) 


73.0(1.5) 


2008 Jun 14 




-0.0209(13) 


+0.0140(13) 


0.0251(13) 


73.1(1.5) 


Overall 




-0.02476(61) 


+0.01706(71) 


0.03007(64) 


72.72(64) 


2008 Jun 08 


Cygnus X-1 


+0.6003(51) 


-4.8315(48) 


4.8687(48) 


138.541(30) 


2008 Jun 09 




+0.7721(34) 


-5.0076(44) 


5.0668(43) 


139.382(20) 


2008 Jun 10 




+0.5467(40) 


-4.9797(29) 


5.0096(29) 


138.132(23) 


2008 Jun 11 




+0.6299(40) 


-4.9701(13) 


5.0099(43) 


138.612(23) 


2008 Jun 12 




+0.6590(40) 


-5.0088(32) 


5.0519(33) 


138.747(23) 


2008 Jun 13 




+0.5537(43) 


-4.9368(46) 


4.9677(46) 


138.200(25) 


2008 Jun 14 




+0.6912(39) 


-4.9267(36) 


4.9749(36) 


138.993(22) 


Overall 




+0.6360(51) 


-4.9579(34) 


4.9985(35) 


138.655(29) 



Note. — TP is measured by APDs (wavelength range of about 400 to 850 nm), and 
remaining targets are observed by PMTs (wavelength range of about 400 to 675 nm). TP 
is subtracted for these targets. Uncertainty in position angle is purely statistical, as position 
angle is not calibrated absolutely. 
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To claim variability from this test, I require a variability confidence level of > 99%. 
Therefore, rejection of the null hypothesis with significance a^s < 0.01 is required to claim 
variability of the source. This method is preferred over a analysis with a constant polariza- 
tion model, because no assumptions are made about intra-night measurement distributions 
and number of measurements (Clarke & Naghizadeh-Khouei 1994). However, K-S tests of 
my data show that such measurement distribution is consistent with a Gaussian nature for 
all targets and for all nights. That is, neither systematic effects nor intra-night variability 
are significant for any targets. Plotted in Figure [3] are the oks values for all pairs of nights 
from each star. The likelihood of variability decreases from Figures [3^ through [3^. 



4.2. HD 189733 (V452 Vul) 

From Figure it can be seen that significant variability of the system is not observed. 
In Figure m I show nightly mean polarization of HD 189733 observed with POLISH and com- 
pared with the data from BOS. It should be noted that the bandpasses are different between 
BOS {B band) and this work (400 to 675 nm), but the discrepancy in Figure H] cannot simply 
be due to this. A analysis shows that the model reported by BOS to fit their observations 
fails to accurately reproduce my observations to the > 99.99% confidence level (x^ = 99.9, 
z/ = 12). Figure [5] isolates the POLISH data, and Figure [6] illustrates the non- variable con- 
trol system HD 175541 for comparison. To find an upper limit to the polarization amplitude 
of HD lS9733b, I perform Monte Carlo simulations for expected exoplanetary polarization 
using Equations [1] through 4. I set i = S5?5S (Torres et al. 200S analysis of Knutson et 
al. 2007), < AP < 3 X 10"^ and 0° < n < 1S0°. There exists a reduced minimum 
of x^/'^ = 1-27 with a = 0.229 probability of a successful fit for AP = 4.0 x 10^^ and 
Q = 169° (Figures [5] and [7]). This amplitude is larger than the maximum AP ^ 2 x 10^^ 
allowed for p = 2/3 (section 2), indicating that it cannot be due to the planet. Indeed, a 
constant fit to the data in Figure [5] produces a similar reduced x^ of 1.79, reinforcing the 
results of the K-S test that my observations of HD 1S9733 are consistent with noise. 

Upper limits to the polarization amplitude can be determined from the cumulative dis- 
tribution function of function of AP (Figure [7|d), which is shown in Figure [SI The 
6S%, 95%, and 99% upper limits to exoplanetary polarimetric modulation are AP < (4.S, 
6.S, and 7.9) x 10^^, respectively. Such large upper limits exist because I was unable to ob- 
serve the system where exoplanetary polarization is maximized, which is near both quarter 
phases. Indeed, the upper limits are based significantly on the single night of observations 
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Fig. 3. — Polarimetric variability of observed targets: (a) Cygnus X-1, (b) HD 189733, (c) 
HD 149026, (d) HD 175541, and (e) HR 8974. Stars with most ors values less than 0.01 
(blue regions) are considered variable. It can be seen that only (a) Cygnus X-1, the variable 
control source, is significantly variable. 




Fig. 4. — Comparison of observed polarization of the HD 189733 transiting hot Jupiter 
system. Solid circles are data from this work, open circles are data from B08, and the 
curve represents the exoplanet model of B08. Mean polarizations from each data set have 
been independently subtracted to show residual polarimetric variability. Phase indicates 
mid-transit in order to be consistent with B08. 
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Fig. 5. — Observed polarization of the HD 189733 transiting hot Jupiter system with best 
fit exoplanet model. Telescope polarization and mean polarization, both effects being at the 
level of three parts in 10^, have been subtracted. Phase indicates mid-transit (inferior 
conjunction of the exoplanet) in order to be consistent with BOS. 
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Fig. 6. — Observed polarization of the long period HD 175541 exoplanet system. Telescope 
polarization as well as residual mean polarization of about 0.1% have been subtracted. 
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at phase 0.65, which is near last quarter (Figure [5]). Since stellar activity inhibits albedo 
constraint from photometry (Rowe et al. 2008a), future observations of this system with 
POLISH are necessary in order to provide an albedo constraint. However, my observations 
have ruled out the AP 2 x 10~^ amplitude of modulation reported by BOS with > 99.99% 
confidence. 

A potential cause of polarimetric variability in hot Jupiter host stars is starspot activity. 
Photometric observations with the MOST satellite suggest the existence of starspots on the 
short period r Boo that follow the rotation period of the star (Walker et al. 2008). There 
is also some evidence that Ca II H and K emission from the short period HD 179949 may 
follow the stellar rotation period (Shkolnik et al. 2005, 2008). HD 187933 itself is known to 
be active, with up to 1% of its surface covered in spots at any time (Hebrard & Lecavelier 
des Etangs 2006; Croll et al. 2007; Pont et al. 2007; Winn et al. 2007; Moutou et al. 2008). 
These spots appear to rotate with the roughly 11.8 day stellar rotation period (Henry & 
Winn 2008; Croll et al. 2008). The run-averaged polarization of HD 189733, observed by 
both B08 as well as POLISH, is an order of magnitude larger than those seen in stars at 
comparable distances (Hough et al. 2006; Lucas et al. 2009). This enhanced polarization 
may be due to starspots. 

Since optical depth of a stellar atmosphere decreases from center to limb, and since the 
scattering angle increases to 90° at the limb, stellar polarization is dominated by the limb 
(Carciofi & Magalhaes 2005, hereafter CM05). Spatially resolved polarimetry is possible for 
the Sun, and this effect has been observed (Faurobert et al. 2001; Faurobert Sz Arnaud 
2003). Just after a starspot appears above the limb and just before it disappears below the 
limb, the greatest asymmetry in limb polarization occurs. Indeed, CM05 model the polarized 
signal from an exoplanetary transit to be strongest at ingress and egress, when part of the 
limb is occulted. The radial polarization profile of most stars is unknown, but the models 
of CM05 indicate that an exoplanet-sized starspot on HD 189733 may generate a polari- 
metric amplitude of 10~^ to 10~^ when near the stellar limb. The unknown stellar latitude 
and stability of these spots must generate complex variations in stellar polarimetry. While 
starspots on HD 189733 appear to be tied to the stellar rotation period, their effects may 
not average out even during long observing campaigns. Unfortunately, B08 do not discuss 
the possibihty of starspots causing their observed modulation. 



0.2 0.4 0.6 0.8 1 

x10-^ 

Fig. 7. — Distribution of (a) x^/ v and (b) as a function of AP for Monte Carlo simula- 
tions. 




Fig. 8. — Cumulative distribution function for Monte Carlo simulations with polarization 
amplitude AP. Confidence intervals and 68%, 95%, and 99% upper limits on AP are shown. 
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4.3. HDE 226868/Cygnus X-1 (V1357 Cyg) 

I verify the exceedingly large significance of polarimetric variability of this system in 
both Stokes parameters (Figure [3^) as well as its roughly phase- locked amplitude of variabil- 
ity of AP pa 0.1% (Table H]). Further analysis of the modulation of this system, including 
inclination estimation, will be left for a forthcoming paper. 



4.4. HD 149026 

A short period, transiting exoplanet orbits this weakly polarized star (Sato et al. 2005). 
Significant variability of this system is not observed to a precision of APobs = 1-6 x 10~^ 
(Figure [3t and Table [3]). Unfortunately, my lack of observations near quarter phase inhibits 
constraint of the exoplanetary albedo. 



4.5. HD 175541 

This weakly polarized star harbors a long-period exoplanet (Johnson et al. 2007). Since 
the expected polarimetric signal due to the exoplanet is APgxp ^ 2 x 10^^ over an entire 
orbit, any observed polarimetric variability from the system cannot be due to the exoplanet. 
This conclusion is bolstered by the lack of significant variability of the system: I observe no 
significant variability above the level of APobs = 2.6 x 10^^ (Figure [3li and Table [3]). The 
trend seen in Stokes Q lies at low significance (Figures [3li and [6]). 



4.6. HR 8974 (7 Cep, HD 222404) 

Hatzes et al. (2003) discovered an exoplanetary companion to the primary component 
of this unpolarized binary system. I expect the amplitude of the exoplanetary polarimetric 
signal over an entire orbit to be APcxp ^ 5 x 10~^ and consequently undetectable. Indeed, 
examination of Figure [3^ and Table [3] shows that this system is not significantly variable to 
a precision of APobs = 9.0 x 10~^. 
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5. Conclusion 

I do not detect the AP ?a 2 x lO^'^ polarimetric modulation of the hot Jupiter system 
HD 189733 reported by BOS with > 99.99% confidence. Using the high precision polarime- 
ter POLISH, I derive an upper hmit to the polarimetric modulation of the exoplanet of 
AP < 7.9 X 10~^ with 99% confidence. I am unable to constrain the albedo of the exoplanet; 
future observations of the system polarization near quarter phases (phases 0.25 and 0.75) 
are required. 

The amplitude of the signal reported by BOS is at least one order of magnitude larger 
than expected given upper limits to the albedos of three other hot Jupiters (Leigh et al. 
2003a, 2003b; Rowe et al. 200Sb; Lucas et al. 2009). In addition, my significant non- 
detection shows the modulation reported by this group cannot be caused by the exoplanet. 
Starspot activity on HD 1S9733 has the potential to introduce polarimetric modulation with 
an amplitude of one part in 10^ to 10^ (cf. CM05). While starspots on HD 1S9733 appear to 
share the stellar rotation period, their unknown latitudinal distribution and timescalc of sta- 
bility suggest that their effects may not simply average out over long observing campaigns. 
However, my data are acceptably explained by a constant polarization model, so polarimetric 
variabihty of the HD 1S9733 system is not statistically significant to the APobs — 2.1 x 10~^ 
level. 

My observations of the transiting hot Saturn HD 149026b encompass neither first nor 
last quarter phases, so I cannot constrain the geometric albedo of this exoplanet. I have 
observed no significant, polarimetric variability of the short period exoplanets HD lS9733b 
and HD 149026b nor of the long period exoplanets HD 175541b and HR S974b. However, 
polarimetric variability of the variable control source Cygnus X-1 is detected with high con- 
fidence. Gaussian distribution of nightly measurements shows that no observed target, even 
the famously variable Cygnus X-1, exhibits significant polarimetric variability during a single 
night. 

High precision, polarimetric monitoring of both transiting and non-transiting, short and 
long period exoplanets is a new field that is poised for major discoveries. The first detection 
of exoplanetary optical, scattered light has been made, and it will soon be detected with 
polarimetry. The atmospheres and accurate masses of close-in exoplanets will be probed. 
The scattering properties of exoplanet host star atmospheres will be studied from polarized 
transit observations, which will be valuable to calibrate high precision, photometric obser- 
vations. The atmosphere of the Earth is effectively removed by high precision polarimetry; 
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therefore, such observations allow ground-based observatories to make important contribu- 
tions to cxoplanetary science that complement the ever-increasing stream of discoveries with 
space-based observatories. 
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